arXiv:1503.08789vl [math.AP] 27 Mar 2015 


Invariant Subspace Method and Fractional 
Modified Kuramoto-Sivashinsky Equation 


A. Ouhadan^ and E. H. El Kinani^ ^ 

^ Centre Regional des Metiers de I’Education et de la Formation, Meknes, 

BP 255, Morocco. 

^A.A Group, Mathematical Department Moulay Ismail University, Faculty of 
Sciences and Technics Errachidia, BP 509, Morocco. 

1 Universite Moulay Ismail Ecole Nationale Superieure des Arts et Metiers 
(ENSAM), Marjane 2, B.P. 15290, Meknes, Maroc. 


Abstract 

In this paper, the invariant subspace method is applied to the time 
fractional modified Kuramoto-Sivashinsky partial differential equation. 
The obtained reduced system of nonlinear ordinary fractional equations 
is solved by the Laplace transform method and with using of some useful 
properties of Mittag-LefHer function. Then, some exact solutions of the 
time fractional nonlinear studied equation are found. 
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1 Introduction 

In the last decade, fractional calculus attracted a great interest of many re¬ 
searchers. The idea of fractional order derivative was started with half-order 
derivative as discussed in the literature by Leibniz and L’Hdpital. Next, it was 
extended to an arbitrary order derivative by Liouville, Riemann, Grfinwald, 
Letnikov, Gaputo etc. In addition, different approaches to define fractional 
derivatives are known [Bnisi SI- The study of fractional differential equations 
becomes of great interest, since for their widely applications including fluid 
flow, dynamical processes in self-similar and porous structures, electromagnetic 
waves, probability and statistics, viscoelasticity, signal processing, and so on 

[mail]. 

The construction of particular exact solutions of fractional differential equa¬ 
tions is not an easy task and it remains a relevant problem. This is the reason 
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why a powerful methods for solving those fractional equations were recently de¬ 
veloped in the literature, including Adomian decomposition method [^, first in¬ 
tegral method [7] , homotopy perturbation method [S] , Lie group theory method 
0 m El m [s] and so on. Most recently, according to invariance princi¬ 
ples, the invariant subspace method developed by V.A. Galaktionov and S.R. 
Svirshchevski m to study partial differential equations was extended by R.K. 
Gazizov and A.A. Kasatkin m to construct some particular exact solutions for 
time fractional differential equations. 

The invariant subspace method used in the present paper, yields us with an 
exact solutions of the time fractional modified Kuramoto-Sivashinsky equation 
in terms of the well known Mittag-Leffler function. In the paper [15], the in¬ 
variant subspace method and Lie group analysis are joined to solve the reduced 
fractional ordinary differential system and the original studied equation. In our 
case, resolution of the reduced system is done by the Laplace transform method 
and by using of some remarkable properties of the well known Mittag-Leffler 

function [isiEiig- 

This paper is organized as follows: In section 2, we recall some main results 
of fractional derivatives and integrals. Section 3, is devoted to describe the 
invariance subspace method. While in section 4, we use the described method 
to construct exact solution admitted by the time fractional modified Kuramoto- 
Sivashinsky equation. Finally, a conclusion is given. 


2 SOME BASIC RESULTS ON FRACTIONAL 
CALCULUS 


This section is devoted to recall briefly some definitions and basic results on frac¬ 
tional calculus. For more details and proofs of the results, we refer to Hiiiiii]. 

The Riemann-Liouville fractional integral is defined by: 



( 2 . 1 ) 


where 7 S R+, and 



( 2 . 2 ) 


is the Euler Gamma function. 

By definition J°/(t) = fit) and it satisfies the property J^Jffit) = f{t). 

First recall that there are various contributions [Diiaii] to define fractional 
derivatives. In this paper, we adopt the fractional derivative in the sense of 
Gaputo [DEIEIII]. The Gaputo definition is used not only because it makes 
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easy the consideration of initial conditions but also because the derivative of 
a constant is equal to zero. In what follows, we recall some important results 
and properties of fractional derivatives and integrals. For more details see for 
example [3]. Before going on, let us denote by AC"([0, t]), n G N the class of 
functions f{x) which are continuously differentiable in [0, t] up to order (n — 1) 
and with/("-!) e AC([0,t]). 


Theorem 1 Letn — 1 < a <n, withn G N. Iff{x) G AC^dOjt]), then the Ca- 
puto fractional derivative exists almost everywhere on [0,t] and it is represented 
in the form: 


Dffit) = Jr"D^f(f) = / [\t-Tr-^-^f^-HT)dT, a^n. (2.3) 

r(n - a) Jo 

The Caputo derivative (12.31) and the Riemann-Liouville integral (I2.ip satisfv the 
following properties [4]: 


JtD?fit) 


= fit), a > 0, 


n —1 






(2.4) 

(2.5) 


Jf‘t'^ 


Dft^ = 


a>0, 7>-l, t>0, (2.6) 

r(7 + 0 + 1) 
r(7 +1) 

—a > 0 ,7 e] - 1,0[u]0, +oo[, t > 0. (2.7) 
i (7 — o + I j 


Here, it is important to mention that the studied equation is a time fractional 
partial differential equation of order 0 < a < 1, so the integer n appearing in 
the relation p.3p is equal to one. Consequently, the formulae (j2.5l) becomes: 


J^D^fit)=fit)-fiO), t>0. (2.8) 


3 DESCRIPTION OF THE INVARIANT SUB¬ 
SPACE METHOD 

The aim of this section is to collect and to present some necessary and essentials 
results from invariant subspace theory. The invariant subspace method m 
was firstly used to construct particular exact solutions of evolutionary partial 
differential equations of the form: 

Ou 

— = F{u,uix,U 2 x,---,Ukx), fc G N, (3.1) 

at 
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where u = u{t, x), Uix = is the i-th order derivative of u with respect to the 
space variable x and F is a nonlinear differential operator. 

Recently, Gazizov and Kasatkin |15] showed that the invariant subspace 
method can be applied also to equations with time fractional derivative: 

In fact, consider the time fractional partial differential equation of the form: 

D'^u{t,x) = F[u], (3.2) 

where F[u\ = F(u, uix, U 2 x, ■ ■ ■, Ukx) and Df is the time fractional derivative in 
the sense of Caputo. The invariant subspace method is based on the following 
basic definitions and results [HIE]. 

Let fi{x),..., fn{x) be an n linearly independent functions and Wn is the 
n-dimensional linear space namely Wn = {fi{x),..., fn{x)). Wn is said to be 
invariant under the given operator F[u\ if K[m] G Wn whenever u G Wn- 

Proposition 3.1 Let Wn be an invariant subspace ofF[u]. A function u{t, x) = 
Sr=i ® solution of equation \3.2i) if and only if the expansion co¬ 

efficients Ui{f) satisfy the following system of fractional ordinary differential 
equations: 

{ Dfui = Fi(ui,...,u„), 

DfU2 = F2(ui,...,u„), 

Ft I Un} — Fn{ui^ ■ • ■ ; 
where Fi,... ,Fn are given by: 

F{cifi{x) H- \-Cnfnix)) = Ki(ci,..., c„)/i(a:) H-h ^^(ci,..., c„)/„(a:). 

(3.3) 

Remark 1 A crucial question in the theory of invariant subspace method was 
how to get the corresponding invariant subspace of a given differential operator. 
This question is solved by the following proposition and for more details see m- 

Proposition 3.2 Let /i(a;),..., fn{x) form the fundamental set of solutions of 
a linear n-th order ordinary differential equation 

L[u] = + ai(x)y("“^) H-h a„_i(x)y' + an{x)y = 0, (3.4) 

and F[y] = F{x, y, y',..., y^'"'>) a given differential operator of order k < n — 1, 
then the subspace Wn = {fi{x),..., fn(x)) is invariant with respect to F if and 
only if: 

L[F[y]] = 0, (3.5) 

whenever y satisfies equation i3.4\)- 
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Remark 2 Condition of invariance appearing in the above proposition is the 
invariance criterion for equation i3.4\ ) with respect to the Lie-Backlund gener¬ 
ator V = F[y]-^. This criterion shows us how the invariant subspace method 
is related to the techniques used in Lie symmetry analysis, see for more details 


4 EXACT SOLUTION OF THE FRACTIONAL 
MKS EQUAION 


In this section, we use the invariant subspace method to construct some exact 
solutions of the time fractional modified Kuramoto-Sivashinsky equation (mKS) 
which is given by: 

L»“(u) =-U 4 a; - + (1 - A)(ita;)^ + 0<a<l, t>0, (4.1) 

where u = u{t,x) and A g] 0; 1[. In the case a = 1, the (mKS) equation (I4.ip is 
a model for the dynamics of a hyper-cooled melt [531121] ■ A more general class 
of such models was introduced and discussed in [25] . 


Proposition 4.1 For any A g]0; 1[ the nonlinear operator F[u\ given by: 

F[u] = -UAx - U2x + (1 - X){Uxf‘ + KuxxY, (4.2) 

admits W 3 = ( 1 , cosyx, sinyx) with 7 = as an invariant subspace. 

Proof. For any function 

/i(t, x) = Cl-b C 2 cosyx-b C 3 sinyx, (4.3) 


with Ci = Ci{t) arbitrary functions, we get: 

F[h] = —y^C 2 cosyx — y^^Ca sinyx-b y^C 2 cosyx-b y^Cs sinyx 

-b(l — A)(y^C 3 cos^ yx -b y^C| sin^ yx — cosyxsinyx) 

-bA(y'^C 2 cos^ yx -b y'^Cg sin^ yx -b 2 y‘^C 2 C 3 cos yx sin yx) 

= (1 - A)y^(C| -b C 3 ) -b (y^ - y'^)C 2 cosyx -b (y^ - ^'^)Cz sinyx S IF 3 . 

Now, we search an exact solution admitted by the time fractional (mKS) equa¬ 
tion (mn) of the form: 


u{t,x) = Ci+ C 2 cosyx -b C 3 sinyx. 


(4.4) 


Consequently, a function u{t, x) of the form (14.41) is a solution of the time frac¬ 
tional (mKS) equation if the expansion coefficients Ci{t) satisfy the following 
system of ordinary fractional differential equations: 


D?C, = (I - A)y2C| + (1 - A)y2C|, 

D-C 2 = y2(l-y2)C2, 

DfCz = y2(l-y2)C3. 


(4.5) 
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To get a non trivial solution needs to assume the condition C 2 ( 0 )C' 3 ( 0 ) ^ 0 and 
for convenience we suppose ( 72 ( 0 ) = ( 73 ( 0 ) = 1. This last condition will be clear 
when the Laplace transform will be used. We start to construct solution of the 
third equation in the above reduced system of ordinary fractional differential 
equations. We mention that, with the Laplace transform it is frequently possi¬ 
ble to avoid working with equations of different differential orders by translating 
the problem into an easy one. 


Recalling some useful properties of the Laplace transform [T] : 

£{i?“/(t)} = s“/(s)-s“-V(0), 0<a<l, (4.6) 

where 

PCX) 

'2{/(i)} =/(s) = / e-‘'*f{t)dt. (4.7) 

Jo 

By putting 0 = 7^(1 — q^) and applying the Laplace transform on both sides of 
the third equation appearing in the fractional ordinary differential system, we 
obtain: 

s“£{C'3(t)} - s“-V3(0) = 0£{C3(<)} , (4.8) 

it yields: 

£{^^3(i)}(s) = ^3^, (4.9) 

then, with the inverse Laplace transform, it gives: 


Csit) = 

where Ea^p{9t°‘) is the Mittag-Leffler function given by: 

nij-ai 

E^,p{en = EZo 


r(m -I- p) 


(4.10) 


(4.11) 


Not that when /3 = 1, E^^i = E^- 

Two last equations in the fractional ordinary differential system (ESI) are 
the same, hence, 

C2{t) = Cs{t) = E^,^{er). (4.12) 

Substituting the obtained expressions of C 2 and in the first equation of the 
system (ESI), it leads to: 

D^Ci = 2(1 - A)72 {E^^^{er)f . (4.13) 

The Mittag-Leffler function does not satisfy the composition property: 

Ea{x)Ea{y) ^ Ea{x + y), (4.14) 
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but it can be observed that the function [milTllIl]: 


= (4.15) 

1 [ai + 1 ) 

does satisfy the composition property: 

E^{x‘^)E^iy‘^) = E^iix + yr). a > 0, (4.16) 

Using of the above relation (|4.16|) . so the equation (14.1311 becomes: 

= 2(1 - X)y^E^y{e{2tr). (4.17) 

Applying J“ on both sides of equation (|4.17|1 . and using integration of the 
Mittag-LefFler function relation [T] (p. 25), we obtain: 

/ {t - = t^+''~^Eay+„ir]t°‘), (4.18) 

^ 1^1 Jo 

where a > 0, (3 > 0 and > 0 , it leads by taking a = v, rj = 02 “ and /3 = 1 to: 

J“U„,i(0(2t)“) = t“£;„,„+i(0(2t)“). (4.19) 

According to the following relation, it yields: 

J“Z)“C'i(t) = C'i(t)-C'i(0), 0<a<l, (4.20) 

we obtain 

Ci(t) = 2(1 - A)7't“i;a.„+i(0(2t)“) + Ci(0). (4.21) 

We assume ( 71 ( 0 ) = 0. Hence, the obtained solution of fractional ordinary 
differential system (14. 5p yields the following exact solution of the nonlinear time 
fractional modified Kuramoto-Sivashinsky equation gH): 

u(t,x) = 2(1 — A) 7 ^t“i?(j,a+i( 0 ( 2 t)“) + Eay{0t^){cos"fx + sin 7 a:), (4.22) 

where 7 = ^ and 0 = 7 ^ — 7 "*. 

5 SOME PARTICULAR CASES 

In this section, we extract some particular cases, precisely exact solutions cor¬ 
responding to A = -A, with m G N* — {1}. 

Case 1 X — ^ 

This particular value of A leads to 7 = 1 and 0 = 0. Consequently, an exact 
solution of the studied fractional equation 63) is given by: 

u(t,x) = —- -— + cosa;-|-sina;, 0 < a < 1. (5.1) 

1 (a -b 1 ) 
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Case 2 A = —, m > 2 

m ’ 


In this case we obtain that 7 = y/m — 1 and 6 = (m —1)(2 —m). The constructed 
exact solution takes the form: 

= 2^-" ~ {E^,i [(m - 1)(2 - m)2"t"] - 1} 

m (2 — m) 

+ [(to — 1)(2 — rn)t°‘] cos \/m — lx 

+ Ea^i [(to — 1)(2 — m)t°‘] sin\/m. — lx. (5.2) 

Now we look for solutions of nonlinear time fractional equation 611 correspond¬ 
ing to a = 1 and a = ^. 

Subcase 2.1 A = —, m > 2,a = 1. 

m. ■ ■ 


According to the relation 

Ei^{z) = E^{z) = e\ (5.3) 

the corresponding exact solution of equation 63) is obtained to be of the fol¬ 
lowing form: 

’ m[2 — m) I J 

+ (5.4) 

Subcase 2.2 A = —, m > 2, a = i. 

m ’ ’ 2 


According to the relation: 

E^^iiz) = Ei{z) = (l + erf{z)-^ , (5.5) 

where 

erf{z) = -^ f eT''^ dy, (5.6) 

Vtt Jq 

the corresponding exact solution of equation 63 is obtained in this subcase to 
be of the form: 




-f 


\/ 2 (m — 1 ) 
m (2 — m) 

(m-l)( 2 -m)t 


|iili ^2(771 — 1)(2 — ~ l| 

(cos y/m — lx + sin y/m— lx.) . (5.7) 


We end this section by giving some corresponding graphs of some found partic¬ 
ular solutions. 
















Figure 3: Solution with a = ^ and Figure 4: Solution with a = ^ and 
TO = 3 TO = 4 



Figure 5: Solution with a 



1 and 


Figure 6: 



Solution with a = ^ and 
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6 Conclusion 


Here, by using the Laplace transform method and some basic properties of 
the Mittag-Leffler function, we succeed to solve the obtained reduced system of 
ordinary fractional equations. Consequently, the invariant subspace method was 
appropriate to construct some exact solutions of the time fractional nonlinear 
modified Kuramoto-Sivashinsky equation. Finally, we note that the method 
used in this paper can be extended to obtain exact solutions of other nonlinear 
time fractional differential equations. 
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